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The linear and nonlinear dynamical susceptibilities of a 
two level system are calculated as it undergoes a transition 
to a decoherent state. Analogously to the Glover-Tinkham- 
Ferrell sum rule of superconductivity, spectral weight in the 
linear susceptibility is continuously transferred from a finite 
frequency resonance to nearly zero frequency, corresponding 
to a broken symmetry in the thermodynamic limit. For this 
reason, the behavior of the present model (the Mermin Model) 
differs significantly from the spin-boson model. The third 
order nonlinear susceptibility, corresponding to two-photon 
absorption, has an unexpected non-monotonic behavior as 
a function of the environmental coupling, reaching a maxi- 
mum within the decoherent phase of the model. Both linear 
and nonlinear susceptibilities may be expressed in a universal 
form. 

PACS numbers: 03.65. Bz, 42.50.-p, 42.50.Lc, 42.65.-k 



I. INTRODUCTION 

A vast body of work has been devoted to understand- 
ing the transition to decoherence in models of a two-level- 
system (TLS) coupled to an infinite set of environmental 
degrees of freedom Q . The central quantity in this tran- 
sition is the dynamical correlation function, or equiva- 
lently, the linear susceptibility, x"(^), of the system de- 
grees of freedom. To our knowledge, no studies of the 
nonlinear susceptibility have been undertaken. Our in- 
terest in the nonlinear response is twofold: 1) experimen- 
tal probes of Macroscopic Quantum Coherence (MQC) 
might access the nonlinear regime and 2) the nonlinear 
susceptibility appears to bear a different relationship to 
the coherence of the TLS than the linear susceptibility. 
A measure of coherence of the system is the quantity of 
spectral weight within the resonance peak of \" { v ) asso- 
ciated with the TLS — this is simply proportional to the 
probability of absorbing a single photon. Absorption of 
a photon at a precise energy is a property uniquely as- 
sociated with the underlying coherence of the quantum 
mechanical system. As the coupling to the environment 
is increased and the system increasingly localizes in one 
state, the probability of single photon absorption dimin- 
ishes. In this paper we address whether two photon ab- 
sorption is qualitatively similar — is it similarly linked to 
the underlying coherence of the system? We are thus led 



to the study of nonlinear susceptibilities of a TLS as the 
system is tuned through a decoherence transition. 

N. D. Mermin investigated a model obtained from the 
well known spin-boson Hamiltonian by keeping only the 
lowest two levels of each harmonic oscillator comprising 
the environmental bath g]. The Hamiltonian therefore 
contains a "system" spin-1/2, represented by the Pauli 
matrices a, coupled to a set of N "environment" spin- 
1/2 degrees of freedom, {sj}. The resulting Hamiltonian 
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(henceforth referred to as the Mermin Model) can be 
solved to demonstrate the close correspondence between 
decoherence and a second order phase transition. The 
transition in the Mermin model also bears some similar- 
ity to the decoherence transition in the Ohmic case of the 
spin-boson Hamiltonian B. Owing to its finite Hilbert 
space, exact dynamical susceptibilities (linear and non- 
linear) can be calculated and followed through the tran- 
sition. The finite dimensionality may also be realistic in 
some physical settings and exhibits the expected correc- 
tions to an "infinite environment" thermodynamic phase 
transition. In this paper, we will present such calcula- 
tions and a careful study of the Mermin Model at the 
linear level, demonstrating an analogy to the supercon- 
ducting transition in the "dirty" limit. We will then go 
on to the calculation and interpretation of the third or- 
der, nonlinear susceptibility corresponding to two photon 
absorption. 

To briefly summarize our findings: The transition to 
decoherence is induced by strong coupling to environ- 
mental degrees of freedom (large A) or a small system 
energy scale (small t). Below a critical coupling, A c , all 
spectral weight of the dynamical susceptibility 



X» = Im 



dt(0\[a(t),a(0)}\0)e(ty 



(2) 



lies in the principal resonance of the TLS at an energy 
t. When the coupling is increased above the critical cou- 
pling, a new exponentially small energy scale 0(te~ N / 2 ) 
emerges, associated with a broken symmetry, (a) 0, 
in the thermodynamic limit N — > oo. This feature 
corresponds to tunneling modified by a Franck-Condon 
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type overlap factor. As the coupling is increased, spec- 
tral weight is continuously shifted to the "near-zero" fre- 
quency channel. The weight of the delta- function, <5(0 + ), 
is simply the order parameter, (0|er|0}| 2 (or, more ex- 
actly, |(l|cr|0)| 2 , for large but finite N). However, the dy- 
namical susceptibility obeys a sum rule implying that an 
incompletely formed broken symmetry state leaves some 
spectral weight at the position of the principal resonance, 
set by t. 

It is in this respect that the Mermin model differs from 
the ohmic case of the spin-boson Hamiltonian (SBH). 
In the SBH, the Toulouse limit (at a — 1/2), corre- 
sponding to complete inelastic broadening of the reso- 
nance, is a precursor to localization of the system spin 
(at a = 1) @,||. In the Mermin model, the quantum res- 
onance of the system spin — although damped — remains 
intact through the decoherence transition. This behavior 
is analogous to the Glover-Tinkham-Ferrell || sum rule 
for the dynamical conductivity in superconductors. Spec- 
tral weight falling in the gap region, proportional to the 
superfluid density (or, equivalently, the order parameter) 
is redistributed to a delta-function at zero frequency. In 
the case of the transition to decoherence, finite but super- 
critical coupling, A > A c , is analogous to the transition 
to superconductivity in the dirty limit where the super- 
conductor can absorb energy above the gap. 

In contrast, the nonlinear response has a more compli- 
cated behavior. Two photon absorption (TPA) requires 
the presence of ancillary levels as intermediate states and 
is therefore enhanced when there is some degree of "in- 
elastic" broadening of the primary resonance. At the two 
extremes, TPA is zero when k > k c and the environment 
is decoupled from the system — but it is also zero when 
k — ► and the system is completely decoherent. 

We now turn to a review of the analytical results ob- 
tained from the Mermin model. 



As Mermin points out, the advantage of the Hamilto- 
nian (^) is that in the limit N — > oo, the environment 
spins may be replaced in the Hamiltonian by the x and 
z components of a classical spin angular momentum: 
jjPzSx — * \Xa z sm9 and ^S z — > — \uocos9. The re- 
sulting Hamiltonian 

H = \^a x + \i z sin 9 - ^ cos 9 (4) 

may be diagonalized and its ground state eigenvalue, 
Eq(0), given by 



E (8) = -- ] Jt2 + - S m 2 9--cos9 (5) 
minimized with respect to 9. 
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II. THE MERMIN MODEL 



Figure 1 



In the absence of coupling to the environment, the 
two level system in (Q) possesses a ground state in a co- 
herent superposition of spin up and spin down (denoted 
|+) — |— )) and will display a sharp resonance in x"i v ) at 
v = t. When environmental coupling is included, there 
are two effects commonly associated with decoherence. 
First, the coherence features should shift to lower fre- 
quency (Franck-Condon overlap) and secondly, the peak 
should broaden from "inelastic" energy exchange with 
the environment. These two effects are nicely demon- 
strated in the Mermin Model. In (|l|), the environment 
spins may be summed to one big O(N) spin, S, and the 
Hamiltonian becomes 



H = — fj T H <t z S t H S z 

2 2N 2N 
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FIG. 1. The classical (N — > oo) minimum energy configu- 
rations of the environment (solid black arrows) in the strongly 



coupled phase, k < k c . (k = 2ujt/\ 
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The critical behavior of the model may now be seen 
by examining the ground state energy Eq, as function 
of t. Eq bifurcates at a finite value of t going from a 
singlet, non-degenerate root for t > X 2 /2ui to a doubly 
degenerate set for t < A 2 /2a;. 







sin 9q = ±1 



'l -4i 2 w 2 /A 4 
l+w 2 /A 2 



t > \ 2 /2lo 
t < X 2 /2u> 



(6) 



In the former case, the environment is decoupled from the 
system and always points along the z-axis to minimize its 
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"Zeeman" energy corresponding to the last term of 
In the latter case, the environment and the system spins 
are frozen in two distinct orientations with degenerate 
energies. The ground state wave function is still in a 
superposition a\+) — 0\—) but it is no longer an evenly 
weighted one (a = /3); rather, the two roots correspond 
to the system predominantly in |+) or |— ). 

However, for finite (even very large) N, the sys- 
tem must be in an evenly weighted superposition — 
independent of the value of t. For t — > 0, it is possi- 
ble to estimate the tunnel splitting between the ground 
and first excited states. Following we consider the 
adiabatic environment (or small t) limit in which the en- 
vironment instantaneously adjusts to the localization of 
the system in a given state and thus "points" at the angle 
6*o, as shown in Fig. 1. For the system to jump to the 
other state, the environment must be caught in a fluctu- 
ation which points along the other direction, —9q. The 
amplitude for N spin-l/2's to be found at an angle 29q 
away is cos^ 9q . The new tunnel splitting is thus reduced 
by a Franck-Condon type overlap factor: 



t* =t 
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, a] R y {29 )\9 ) =;cos iV e Q =t 
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where R v (29q) is the rotation operator. Note that this 
is still a purely conservative process; the resonance is 
perfectly sharp although at an exponentially reduced en- 
ergy scale. We will refer to this resonance as the Franck- 
Condon resonance. 

Although, for finite N and t, the symmetry is unbro- 
ken and the ground state is an evenly weighted super- 
position, there must be some symptom of the crisis at 
t ~ \ 2 /2ui when N becomes large. Presumably, this 
symptom should be the inelastic broadening of the res- 
onance to an over damped state. Following along the 
lines of the Fermi Golden Rule (FGR) calculation in [Q , 
we treat the last two terms of (||) as the unperturbed 
Hamiltonian, Hq. Hq is a Zeeman-like Hamiltonian with 
a magnetic field that depends upon the z-component of 
the system spin. Therefore it should be diagonalized with 
a a z dependent rotation about the y-axis: 
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R = Ry{(T z 8()) - 

The transformed Hamiltonian is then 



(9) 



H' Q = R^HqR = ^-S z cos# (l + 4r) (10) 

2N LO z 



where tan#o = X/lu has been chosen to eliminate the S x 
term in H' . Now turning to the first term of (||), the 
rotation yields: 

H[ = ^a x R = t -(a+e- la ' e ° s y + a~ e +ia * eoS v) (11) 



Now the system is prepared in a state |+) with the en- 
vironment in its ground state S z = —N/2 (which will 
be denoted | — N/2)). The transition rate is now calcu- 
lated by applying Fermi's Golden Rule, treating H[ as 
the perturbation. 

2 N/2 

T = T £ |(-|Ha- e +-=^|-iV/2)|+)| 2 ( 5( J B / - J B. J ) 

m=-N/2 

(12) 
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Figure 2c 





fore shift from O(i) to a smaller energy scale Ofte-*/ 2 ). 
In the latter case, the system forces the environment 
to dissipate energy and the resonance is broadened by 
0(X 2 /uj) but remains nominally at 0(t). 




FIG. 2. Sequence of states (a) initial, (b) intermediate 
and (c) final, depicting Fermi's Golden Rule applied to 
the Mermin Model. Solid black arrow is the environment; 
state of the system is shown in the upper left. Final en- 
vironment state must make a fluctuation of approximately 
Am = — JV /2cos#o + N/2 to have a reasonable overlap with 
the intermediate state. The system drives the environment to 
dissipate energy. 

As depicted in Figs. 2(a-c), H[ induces a rotation of 
the environment by —9q about the y-axis as the system 
flips from |+) to |— ). Since the environment is now ro- 
tated away from the environment ground state of the final 
configuration ((— 1(— N/2\) as well (Fig. 2c), the domi- 
nant contribution to T comes from an excited (rotated) 
final state. Looking at the geometry of Fig. 2c, the final 
state, (— |(m = —N/2 cos 0o\, contributes most strongly. 
Using the Hamiltonian Hq, the energy difference is com- 
puted: 

uj A 2 

AE = Et — Ei (l-cos6>o)cos6Ul + A 2 /^ 2 ) 

4 8ll> 

(13) 

This calculation suggests an interpretation of the ear- 
lier result (JsJ) for the critical coupling obtained in the 
N — > co limit. When the width of the resonance AE 
becomes comparable to t, the resonance is over damped 
and quantum coherence is destroyed. 

Although both the adiabatic computation and the 
FGR one are small t perturbation theories and involve 
the same environment overlap amplitude, their interpre- 
tations are quite different. In the former, the system in 
the |+) state leaves an "imprint" upon the environment 
which, because it is nearly orthogonal to the imprint left 
by system state |— ), reduces the tunneling amplitude by 
an exponentially small factor. The resonance must there- 



Figure 3 

FIG. 3. Possible "continuous" evolution of spectral weight. 
The principal resonance and the Franck-Condon resonance 
are the same. 

One might speculate upon two possible scenarios in 
which these two perturbativc calculations are blended, as 
suggested in Figs. 3 and 4. In Fig. 3, a single resonance 
broadens and shifts to lower energy, decoherence occur- 
ring when the resonance becomes broadened beyond its 
energy scale. In this scenario, there is no explicit Franck- 
Condon energy scale, just the appearance of increasing 
spectral weight close to zero as the resonance becomes 
over-damped. In the other scenario, demonstrated in Fig. 
4, spectral weight is only transferred to an exponentially 
small energy below some critical coupling; the resonance 
at t, although broadened, is not critically damped and 
remains distinct from the Franck-Condon resonance. To 
put it another way, the systems retains a gap of 0(t). 
Furthermore, even at finite N, there is essentially no 
"Franck-Condon effect" for k > k c (see Fig. 5). It is the 
second scenario that is manifested in the Mermin Model. 



III. CALCULATIONS OF THE LINEAR 
SUSCEPTIBILITY: x (y) 

We consider the Hamiltonian ([!]) with an additional 
coupling to an external field h(t) — hcosvt: 



H ext = -h—cosvt 



(14) 



For calculating nonlinear susceptibilities it will be nec- 
essary to use explicitly real driving fields, so we adopt 
explicitly real notation for the expectation value of the 
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time dependent system spin, at this point. The time de- 
pendence of a z , at linear response, is 



1 



(u z (t)) — hxiy) cos ut — h\' {v)sva.vt 



(15) 



At zero temperature, the response function xi v ) ' IS given 
by the integral on the right hand side of (||). Expanded 
out, 
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80 environment spins. Computations with up to 300 en- 
vironment spins suggest that there is little change be- 
yond 80 spins. An artificial broadening (8 = 0.01) was 
added to make the features more visible. As seen in this 
sequence, n = n c is marked by the appearance of the 
Franck-Condon resonance at an exponentially small en- 
ergy scale (vio ~ 10 -3 ). The spectral weight remaining 
at 0(t) when k > k c clearly exhibits inelastic broaden- 
ing, although the resonance essentially disappears before 
becoming critically damped. 



where j/y = (Ei — Ej)/h and S = + . The matrix ele- 
ments of the system spin cjy = (i\cr z \j) will be referred to 
as the order parameter matrix. (The state \i) now refers 
to the combined system and environment state.) Using 
( |l5| ) , the rate of energy absorption is then 



(17) 



Equations analogous to (jl||), (|l6| ) and ([l7|) for the third 
order nonlinear susceptibility are derived in the Ap- 
pendix. 

The Hamiltonian (Q) was diagonalized and the com- 
plete order parameter matrix, <7y, necessary for the lin- 
ear and nonlinear calculations, was computed. It is nat- 
ural to work with one dimensionless coupling constant, 
k = 2ujt/\ 2 . The critical coupling implied in the N — > oo 
calculation is then k — n r = 1. 
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Figure 5 

FIG. 5. The order parameter, aio, is plotted as a func- 
tion of « and exhibits second order phase transition behavior. 
k was varied by changing t. Comparison of N = 80 and 
N = 300 data show the slope steepening at the critical point. 
When k > k c , there is no weight in the Franck-Condon reso- 
nance other than that attributed to the rounding of the phase 
transition at finite N. 

The counterpart, at finite A, to (a) ^ at infinite N 
is the spectral weight at the Franck-Condon resonance, 
(7xo 7^ 0. Fig. 5 demonstrates that cio behaves as ex- 
pected in a second order phase transition; the slope at 
n = grows with increasing N. The matrix, <zy, and 
consequently, x"{ v )i obey a sum rule: 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.i 



Em 2 = 1 



(18) 



Figure 4 

FIG. 4. Evolution of spectral weight in the Mermin Model. 
x"{ v ) i g shown for four values of k decreasing below the criti- 
cal value, k c . Despite some change in the position of the prin- 
cipal resonance, spectral weight emerges at v ~ abruptly at 
k c and the two features remain distinct. 

Fig. 4 shows the behavior of x"{ u ) as the system en- 
ergy scale t is reduced, corresponding to a range n = 
0.6k c — 1.2k c . These computations were performed for 



An incompletely developed broken symmetry (|ctio| 2 < 1) 
means that spectral weight remains at the principal res- 
onance at 0(t). This behavior is somewhat analogous to 
the superconducting transition and the Glover-Tinkham- 
Ferrell sum rule for the real part of the dynamical con- 
ductivity, c(w): 



(Tn(w)dw 



7ms e 
2m 



<Js((jj)du) 



(19) 



(Subscripts "n" and "s" refer to "normal" and "super- 
conducting".) The sum rule states that spectral weight 
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missing from dynamical conductivity above the gap must 
reappear as superfluid (DC) Drude weight. This behavior 
may be compared with other models and computations 
of environmental spin decoherence S . These models dif- 
fer from the Mermin Model in two respects: 1) typically 
the environmental coupling constant is comparable to the 
system energy scale instead of weak 0(X/N) coupling, as 
in the present case; 2) the coupling between the system 
and environmental spins is isotropic, of the form, a ■ Sj. 
These differences lead to qualitatively different decoher- 
ence effects. 



In the expression above, the energy denominators have 
been abbreviated, 
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it 



(21) 



To calculate (|20|) it is necessary to find all possible sets of 
four couplings that connect the ground state back to the 
ground state. These are referred to as "four link chains" ; 
a few such chains are depicted in Fig. 8. 



IV. NONLINEAR DYNAMICAL 
SUSCEPTIBILITY 

In studying the nonlinear susceptibilities, we are inter- 
ested in power absorption, as this is the relevant experi- 
mental probe for MQC. Consider the response of just a 
single two-level-system (TLS) coupled to ancillary levels: 
Driving the system with a signal at frequency v results 
in third order nonlinear responses at v and 3ia The sig- 
nal at 3v is not absorptive and, in any case, would not 
be detected in any experimental setting where pumping 
and detection take place at (nearly) the same frequency. 
If the nominal frequency scale of the TLS is vq, there 
is an imaginary but emissive response at vq. (The non 
vanishing, odd order, nonlinear corrections to power ab- 
sorption at vq must alternate in sign to produce finite ab- 
sorption at i/q, consistent with the exact Rabi solution.) 
On the other hand, third order two-photon absorption 
(TPA) process at vq/2 would resemble ordinary linear 
absorption of a i/q/2 TLS in that it would be detected as 
an absorptive phase shift when the pumping frequency 
reaches vq/2. The reader is referred to ref. for further 
details about nonlinear processes. 

Fig. 6 shows the two photon absorption spectrum 
along with the linear susceptibility for several values of re. 
Our basic result is the following: as the critical value of 
re is reached and the two-level-system undergoes a tran- 
sition to decoherence, the two photon absorption power, 
^27, begins to increase. It reaches a maximum at approx- 
imately re = re c /2 and then decreases to zero as re goes to 
zero. Fig. 7 shows the behavior of the TPA signal as a 
function of the control parameter re. 

The two photon absorption spectrum was calculated 
from an exact diagonalization of the Mermin model using 
the following expression for TPA (which is derived in the 
Appendix) : 
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Cl 2l = -hW 3] (v) 
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Figure 6 

FIG. 6. Nonlinear susceptibilities for two photon absorp- 
tion compared with the linear susceptibility for two values of 
re. Note the increase in x' 3 ' as re is decreased from 0.8re c to 
0.6re c . 

To understand the non-monotonic behavior of the TPA 
signal, we look at the expression (^0|). Isolating the ab- 
sorptive pole, x^(y) can be factored in the following 
way: 



where a m (z/) is a sum over intermediate states, 

&0n @nrn 



(22) 
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OLmiy) is known in nonlinear optics literature as the 
degenerate hyperpolarizability. We focus our attention 
upon the integrated spectral weight, shown below for 
both linear and nonlinear susceptibilities (fi = 1): 
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x\v)dv 



nonlin 
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(24) 



The limits of integration are meant to restrict the inte- 
gration to the principal resonance and exclude the contri- 
butions from the near degenerate ground states. (These 
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are the states m — 0, 1 excluded in the sum below.) The 
integrated weight is in some sense a "figure of merit" for 
detection of an MQC resonance. A meaningful compari- 
son between linear and nonlinear absorption requires the 
multiplicative factor h 2 in the nonlinear expression. In 
both cases, the integration collapses the absorptive pole 
leaving, in the linear case, 



7 lin = n{ \) 2 E \ a ° n 
and in the nonlinear case, 



*(i) 2 (l-*Si) 



Wilin = 1Th2< ^2 )A % E l a "K>o)| 5 



(25) 
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Figure 7 

FIG. 7. Integrated spectral weight I non \[ n for two photon 
absorption as a function of k. Solid line is the result of theory. 
The applied AC field strength h = 1. 



Since ijj n given in equation (25) depends only upon the 
order parameter ooi , it has a universal behavior in terms 
of k through the decoherence transition (see equations 
(|35|), (pitf)). It turns out to be possible to express I non \[ n 
in a similar fashion within a simple approximation and 
understand the non-monotonic behavior seen in fig. 7. 
It will be shown that I non \[ n depends upon the order 
parameter through: 



nonlin 



^0i) + O(-r) (27) 



and thus reaches a maximum at approximately n = k c /2. 
(In the exact solution, |(0|er z |0}| 2 ~ 1 — n 2 for k c = 1; see 
equation ([35]) .) 

To understand this result, consider the usual condi- 
tions through which a broken symmetry is accessed math- 
ematically: Apply a small static bias field h^c coupled to 
the system through H^ c = —hd c (j z and take the limits in 
the order, 



lim lim 

/idc— >0 N—nx 



l°*|0> 



(28) 



In the decoherent phase [n < k c ), a field /i<j c ~ 0(te~ N / 2 ) 
will localize the system yielding (0|cr|0) ^ 0, but will have 
essentially no effect in the coherent phase (k > n c ). The 
relevant energy levels now look like the right half of Fig. 
8 (the quasi-degenerate ground state levels are no longer 
coupled through a z ). 





Figure 8 



FIG. 8. Two possible four-link-chains. 

The key observation is that the net transition am- 
plitude (proportional to a m {v)) to go from |0) to a 
state \m) through degenerate intermediate states \n), \l), 
etc. is zero unless there is a broken symmetry leaving 
(0\<j z \0) 0. If the intermediate levels are degenerate, 
a m (v) may be written 



E 



1 



V 



VnO 



E 



-((0|a z |0)(0|a z |m) + (m\a z \m)(0\a z \m)) 



(29) 

where the last line is obtained by inserting a complete set 
of states into (0|(t 2 |to} = 0. Since the expectation value 
of a z is 0(1) only in the ground state (in the broken 
symmetry phase), we arrive at 



■|<T z |0)(0|<r z |m) 
v - v n0 



V 



(30) 



Unless (Too is non zero, the transition amplitudes for 
0) — > \m) involving an intermediate state interfere de- 
structively. Finally, the integrated weight is obtained 
from (|2T 
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J nonlin = const x £ \ a ™("™o)\ 2 

m^O 

const ^ 2 

= : ^ X 2—/ a 00 a 0rn&m0 



const 
(y - v n0 f 



X O-QoC 1 - °Oo) 



(31) 



where the sum rule (|l8|) (with j = 0) was used in the last 
line. This is the argument used to establish the result 

We now fill in the details of this argument appropriate 
for the numerical computation performed in this work 
by 1) setting h^c = and 2) including the finite width 
of the resonance. a m (y) is to be evaluated at the fre- 
quency for TPA, v — — this frequency is about 
half of the typical frequencies of the intermediate levels 
v n Q. a is rewritten in such a way as to separate the quasi- 
degenerate ground states from the excited states at the 
energy of the principal resonance at 0(t) (see Fig. 9): 



t (^ m0 /2) = 



+ £ 

n#0,l 



(32) 



^io = t* ~ 0(te~ N / 2 ) is the energy scale of the Franck- 
Condon resonance and will be set to zero. The difference 
between the resonant frequency, u m o, and the intermedi- 
ate state level, is n o, which appears in the second term of 
( |32| ) is bounded by the width of the primary resonance 
and is therefore typically much smaller than the primary 
resonance itself. Denoting a maximum value for that dif- 
ference by e (that is, is n o — I'mo + e )j the second term of 
(E52I) becomes: 



&0n &nri 
^mo/2 



y 



(33) 



where vq is the nominal energy scale of the resonance. 
Following the previous calculation for ha c ^ 0, the hy- 
perpolarizibility a is now 



40-oiCTim e . 

a m {vo) = h 0( — ) 



(34) 



Substituting a back into equation (|26|) and using the sum 
rule ( |l8| ) (with j = 1) establishes the result (p7|). 

To express and ^ non i m in terms of the control 
parameter k, we appeal to the exact solution in the N — > 
00 limit. Specifically, the matrix element ooi is related to 
the expectation value of a z in the either of the degenerate 
ground states of (0) (for k < k c ) by: 



lim |cr i| = |(cr z ) | = 

N— »oo 



1-K 2 



1 + Kn 2 \ 1 + -« 



1 — K 2 



Therefore and ^ non i m may be compactly written in 
terms of n: 



(36) 



_ 2nh 2 1 4 1 - k 2 1 - k 2 

'nonlin = -jfty TT4^ (1 " 1 + K* 2 ' (3?) 



p states 



even 




/ odd 



Figure 9 



FIG. 9. Depiction of the quasi-degenerate ground state 
|0), |1) coupling to p excited states. The splitting between 
|0), |1) defining the Franck-Condon resonance is 
The set of closely spaced excited states represents the inelas- 
tic broadening of the principal resonance. 

Returning to fig. 7, the non-monotonic behavior of 
^nonlin as a function of K follows from its dependence 
upon the combination: (er 2 ) 2 (l — (a z ) 2 ). A plot of ^ n0 nlin 
following from the analytic expression ([37]) is shown for 
comparison. In the numerical computations, the max- 
imum intensity appears at approximately k = 1/2 in 
agreement with theory. 



V. CONCLUSION 

In this work we have attempted to make more explicit 
the relationship between the decoherence transition in 
the Mermin Model and a typical thermodynamic phase 
transition. At the onset of decoherence, spectral weight 
is transferred from above the "gap," set by the system 
energy scale t, to a delta function at v — 0. This behav- 
ior is in contrast to the SBH where symmetry breaking is 
prefaced, first by the complete disappearance of the in- 
elastic peak in x"( v ) & t finite frequency (when a = 1/2), 
and then finally by the softening of renormalized system 
energy scale to zero (when a = 1). 
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Also, we have studied nonlinear absorption in a deco- 
herent TLS. Single photon absorption is a measure of 
coherence of the system because it samples the spec- 
tral weight remaining at the principal resonance — or, 
the complement of the order parameter according to the 
sum rule. On the other hand, two-photon-absorption 
(TPA) samples the product of the order parameter and 
its complement and thus vanishes in the fully coherent or 
fully decoherent phase. TPA reaches maximum intensity 
somewhere between these two extremes. 

Observing TPA may be possible in systems demon- 
strating quantum coherence such as the recently studied 
Fes molecular nanomagnet. In principal, measure- 
ment of the ratio of I non \[ n to /j m would provide a di- 
rect measure of the order parameter (0|<r z |0), the local- 
ized "weight" of the system. If we consider a sample in 
which the system spins are individually partially local- 
ized but uncorrelated with each other (like a spin glass), 
TPA may be a useful probe for this type of order. 
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APPENDIX A: DERIVATION OF THE 
NONLINEAR DYNAMICAL SUSCEPTIBILITY 

We begin by expressing the full time dependence of the 
expectation value of the desired operator in the presence 
of the driving field. 

(cr,(t)> - (0\U\t, -oo)a,(t)U{t, -oo)|0) (Al) 



where 



U(t, -co) = Te 



-i p H(t')dt' 



1 -i 



dtiH(tx) 



dt x / dt 2 H(tx)H{t 2 ) 



+i 



dti 



h 

1 dt 2 I 2 dtzH(ti)H(t 2 )H(t$) 



(A2) 



The time dependent operators appearing in these expres- 
sions are interaction picture operators, treating the full 
Hamiltonian (|l|) as the bare Hamiltonian and the exter- 
nal coupling Hamiltonian , H(t) — —^cr z (t) cos ft, as the 
perturbation. The S'-matrix is expanded to the appro- 
priate order. Grouping the third order terms in (Al) 
together, we arrive at: 

(0Ki)/s|0> + (0|4ff(t)|0> - 

dt x / dt 2 / dt 3 (0|cr(t)cr(ti)<r(taMt3)|0) 



xcosi^i cos vt 2 cos vt-jfi(t — ti)8(t\ — t 2 )6(t 2 — t^) 

+ c.c. (A3) 

and 



(0|/ 1 t a(t)/ 2 |0> + (0|7 2 t a(t)I 1 |0) = 
+ih 3 



dti 



dto 



dt 3 (0|cr(*i)cr(t)<r(t2)(7(f3)|0) 



xcosi^i cos vt 2 cos vt?fi(t — t\)Q(t\ — t 2 )6{t 2 — t^) 

+ c.c. (A4) 

We use the integral representation of the step function, 



i r 



(ko- 



to—id 



(5 = + ) and collapse the time integrations. From the 
various products of delta functions, terms are grouped 
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according to whether they oscillate at ±3^ (third har- 
monic generation) or ±i/. From the latter terms, absorp- 
tive contributions come from the part of the signal that 
is 7r/2 out of phase with the driving term, — \a z cosvt. 

Terms that modify absorption at the original linear 
resonance, x" ~ P n 'o{ v ) can be grouped as follows: 

nml 

x[i^ l (-i/)i^ m (-2i/)J^ l (-i/) 

+PL(-^)PLmW - W)] (A5) 

The two photon absorption terms can also be grouped as 
follows: 

^27 = T ( 77 ) 4jy TT &0nPnm<7ml & 10 
4 2 ^ ^— ' 

nml 

xP^ m (-2z/)P/ ( l ,)P; ( ! /) (A6) 
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